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Abstract 

In a previous paper we have introduced a class of multiplications of dis- 
tributions in one dimension. Here we furnish different generalizations 
of the original definition and we discuss some applications of these 
procedures to the multiplication of delta functions and to quantum 
field theory. 

Resume Nous avons introduit dans une publication precedente une 
classe de multiplications entre distributions a une dimension. Nous 
donnons ici des generalisations differentes des definitions originelles et 
nous discutons des applications de ces methodes a la multiplication 
des functions delta et a la theorie quantique des champs. 



1 Introduction 



In the past years many attempts have been done to extend the ordinary 
multiphcation between functions to distributions. Of course, as for the ex- 
tensions of (almost) any kind, the procedure is not unique and, in fact, many 
inequivalent proposals are nowadays present in the literature, see [1, 2, 3] 
among the others. 

In this paper we generalize a class of multiplication of distributions in- 
troduced in a previous work by the author, see [4]. The original definition 
was based on two different regularizations of distributions, the analytic and 
the sequential completion methods. In particular, this last procedure makes 
reference to functions in 7){R) which generate the so-called delta-families. A 
delta family is essentially a set of functions which approximate 6{x) in the 
topology of V{R). Sometimes, whenever the applications require it, it may 
appear necessary to use a weaker form of this procedure. Possible weakening 
of the requirements in [4] are part of the containt of this paper. In partic- 
ular, in Section 2 we relax some of the requirements given in [4], so that, 
in principle, more distributions can be multiplied between themselves, while 
in Section 3 we give two inequivalent multiplications among more than two 
distributions. 

The necessity for extending the multiplication to more than two distri- 
butions follows directly from physical examples: this is what we need to 
do whenever we try to regularize three or four-points Green's functions in 
a given quantum field model. In fact, it has been recognized since Wight- 
mann's work, [5], that the field operators are not operator- valued functions, 
but rather distributions defined on a certain domain dense in a Hilbert space. 
We know also that a field theory is often defined via a lagrangian density, C, 
which depends on the products of such fields considered at coincident points. 
This is, of course, an operation which has no rigorous mathematical meaning 
in this naive form. One of the most famous consequences of this procedure 
is that certain Feynman diagrams diverge, and that analogous divergences 
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are observed also in many matrix elements of the dynamical variables. Many 
attempts have been made in the past decades to give a rigorous meaning 
to quantum field theory (QFT). This has generated essentially two different 
approaches: constructive QFT, as proposed first by Wightmann, where the 
Wightmann functions, that is some 'matrix elements' of the fields, are the rel- 
evant dynamical variables, and the algebraic QFT, developped by Haag and 
Kastler in the sixties, which we will not consider here. In [5] it is widely dis- 
cussed a sort of regularization procedure for the Wightmann functions: they 
can be recovered as the boundary values of some holomorphic functions. Nev- 
ertheless, even using this regularization, the problem of the divergences in a 
perturbative QFT still exists and the solution proposed are, in our opinion, 
unsatisfactory. 

A possible way out is to compute these Feynman graphs using the regu- 
larizations of the fields, instead of the fields themselves, and to remove the 
regularization only at the end. This is the path we will follow in Section 
5. The results we will obtain, however, show once more the difficulty of the 
problem: in particular we will see that our regularization procedure, as it is, 
is not powerful enough to avoid the appearance of the divergences in a free 
QFT in 1+1 dimensions. 

An idea close to our is also behind Colombeau's book, [6], where a sys- 
tematic approach to quantum field theory (QFT) is proposed. The lack of 
uniqueness in the regularization procedure makes Colombeau's work not res- 
olutive. For instance, other approaches, more along the lines of this paper, 
can be found in [2, 7]. 

The paper is divided as follows: 

in the next Section we start recalling the definition of the multiplication 
given in [4]. We take also the opportunity for briefiy discussing some new 
result. Then we 'relax' this definition of the multiplication to better deal 
with physical models involving distributions of <S', like in QFT; 

in Section 3 we discuss two possible generalizations of the theory discussed 
in Section 2 to more than two distributions; 
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in Section 4 we show many examples involving delta functions of both 
the generalizations proposed. We also discuss some physical examples; 

in Section 5 we apply our procedure to an easy quantum field model, the 
free 1 + 1 Klein-Gordon theory; 

in the last Section, finally, we comment the results and state a plane for 
the future. 



2 Definition of the multiplication 

In this Section we briefly recall, for readers' convenience, the basic definitions 
and results of the multiphcation introduced in [4]. We will shghtly modify 
this definition in a second time, so to build up a framework which is more 
suitable for apphcations to QFT. 

Let V be the subspace of all the functions in C°° with arbitrary support, 
8, with the following properties: 

i) 0(x) |x| < /co for |a;| — > oo, 

ii) \x\ < kn for a; — > oo, 

where ko, ki, ... are constants. The convergence is defined as in S. 
Let V be the dual space of V. For distributions in this space it has been 
shown in [2] that the function 

T\z) = —T-(x-z)-^ (2.1) 

exists and is holomorphic in z in the whole 2;-plane minus the support of T. 
The function 

Tred{x, e) = T°(x + ie) - T°(x - ie) (2.2) 

is further weakly convergent to the distribution T when e goes to zero, [2]. 
Also, if T{x) is a continuous function with compact support, then Tred{x, e) 
converges uniformly to T{x) on the whole real axis for e — > 0+. 

The other ingredient of the multiplication in [4] is the method of the 
sequential completion, which makes reference to the so-called 5-sequences. 
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In [4] we have called (5-sequence a sequence of functions Sn{x) = n(l){nx), 
where (p e I^(R) is a given function with supp C [—1, 1] and / (f){x) dx — 1. 
Then, VT e V'{R), the convolution r„ = T * 5„ is a C°° -function, for any 
fixed n e N. The sequence T„ converges to T in the topology of T>', when 
n — > oo. Moreover, if T{x) is a continuous function with compact support 
then Tn{x) converges uniformly to T{x). 

In [4] we have proceeded in the following way: 
for any couple of distributions T,S e V, Va, /3 > and V^^ E T> we have 
defined the following quantity: 

{S T)("'/^)(*) ^ - / [Si^\x) TUx, -) + Ti^\x) SUx, -)] "^{x) dx 
1 J-oo vr vr 

(2.3) 

where S^^\x) = {S* ^^^^){x), with b^^\x) = n/^^in^x). 

Hence, the two distributions S and T in V' are said to be multipliable if the 

limit of {S ®T)l^'f^^ for n ^ oo exists finite in a weak sense. Finally, we have 

defined 

{S ® T)(«,^)(xl/) ^ Jim (5 ® rt'^\^). (2.4) 

In [4] we have proved, among other things, that this product extend the usual 
product of the functions, in the sense that if T{x) and S{x) are two contin- 
uous functions with compact supports then the product Tjjl^\x) Sred{x, ^) 
converges uniformly to T{x) S{x). As mathematical applications of our def- 
inition we have discussed the possibility of multiplying two (derivatives of) 
delta functions localized at the same point. Here we want to make this in- 
formation complete. In particular we want to extend the multiplication to 
arbitrary derivatives of 6{x), 5^''\ Since this result is a straightforward gen- 
eralization of what has been done in [4], we will not give all the details. If 
we want to define {d^''^ ® 5*-'^ )(«,/?) then we are forced to consider different 



5 



situations depending on the parity of the integers k and I. We get 



a>{k + 1 + 2)13 

i^±^Afc+z+25, a = {k + l + 2)P k,l even 

_{k+l±}^Ak+i+25, a = {k + l + 2)(3 k,l odd 

a = {k + I + 2)/5, k even and / odd, 



(2.5) 



where we have taken, as in [4], 



S^-exp{^}, < 1 
0, > 1. 



(2.6) 



Here F is a normahzation constant, and we have defined, whenever they 
exist, Aj = dx, for integers j. In order to have a finite regularization 

(2.5), we always have to choose a function $(x) with m even and such that 
m > k + I + 1. It is interesting to notice that equation (2.5) imphes, among 
the others, the following equalities: 

{5" 5")(a,f3) = -iS' 5"%,p) = {d® d"'%,p) 
(5'®5')(a,/3) = -(5®^")(a,/3), 

which show that, at least for this particular example, the usual property of 
the derivatives of the distributions are satisfied by the product ®(a,/3). 

As we have already discussed in the Introduction, we are interested in 
applying our proposal of regularization and multiplication to quantum fields, 
which are operators whose matrix elements belong to S'. It is therefore nat- 
ural to generalize a bit the above definition, trying to construct a framework 
more directly related to the physics. In particular, we modify the definition 
of the sequential completion, which is strongly related, in its original version, 
to distributions in V. There exist also technical reasons which suggest to 
relax the definition of the sequential completion. We will comment on this 
point in Section 5. 
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Using the terminology of [8] , we call delta sequence of Dirichelet type a 
family of functions 5k{x), which satisfy the following conditions: 

rA 

i) / Sk{x) — > 1 when /c — > oo for a certain A> 0; 

J -A 

+ / ){5k{x) f{x)dx) = 0; 

c 

in) 3Ci,C2, positive constants : \Sk{t)\ < -rj + C2. 

In particular, it is an easy exercise to prove that for any belonging to 
S, such that ^{t)dt^l, then the family of functions S^^\x) = 
generates a delta sequence of Dirichelet type if /3 > 0. 

We have the following: 

Proposition 1. 

Let T e S'iK) and 6i^'\x) be a delta sequence of Dirichelet functions. 
Then the convolution T^^-* = T * 6ll^^ is a C°°— function, for any fixed n G N. 
The sequences 5^^^ (x) and X^^^ (x) converge respectively to 6 and to T in the 
topology of S', when n — >• 00, for all P > 0. 

Moreover, if T{x) is a Holder-continuous function with compact support 
[a, b], then T^f^^x) converges uniformly to T(x) on every interior subinterval 
of [a, b]. 

Proof 

The C°°— nature of the function Tlf \x) is a well known property of the 
convolutions, which can be found, for instance, in [9]. 

In reference [8] it is proved that any 5— sequence of Dirichelet type weakly 
converges to the 5 function with respect to any function which has a finite 
derivative in the origin. A fortiori therefore will converge to 5 in the 
topology of S'. From this fact it easily follows the convergence of Tjf^ to T 
in S'. 

The last statement is again contained in [8]. □ 
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Remcirks .- (a) Another generalization of the (5-family is also discussed in [8]. 
The family is now called delta family of positive type and, as the name itself 
suggests, its functions must all be not negative. This is, in general, a strong 
requirement which is not necessarily satisfied by the 'generating' function 
we will use in the application to QFT, and this is the reason why we 
have focused our attention to Dirichelet's type functions. Nevertheless, even 
for such a delta family a Proposition like the one above can be stated; minor 
differences are required in the hypotheses but the results, essentially, coincide. 

(b) Any $(t) G V{R) generating a 'standard' delta family also generates 
a delta family of Dirichelet type and, if ^{x) > for all x & R, also a delta 
family of positive type. 

(c) One may wonder why we have introduced so many families of delta 
functions: the reason is that the choice of the function $ cannot be made in 
general a priori by us, but it is often forced by the model which has to be 
regularized. In particular, in the example of QFT we will show that there is 
no reason, in general, for $(x) to have compact support or to be positive. 

(d) Proposition 1 can be used to show that the new multiplication still 
extends the usual multiplication of continuous functions with compact sup- 
port, in the sense that if T{x) and S{x) arc Holder-continuous functions with 
compact support in [a, b] then, Va, (3 > and V^' G 5, then 

/oo 
T{x)S(x)^f(x) dx. 
-oo 

The definition of the multiplication is now, formally, the same as in eqs. 
(2.3) and (2.4). The only difTcrcncc is in the mathematical nature of 

We end this Section giving the extension of the definition of multiplication 
to the case in which the distributions S and T do not commute, even if we will 
not meet with this problem in this paper. In this condition we are forced to 
symmetrize the original definition (2.3), (2.4). Let S and T be two operator 
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valued distributions. Keeping the same notation as before, we define 

1 1 1 

{S ® r)(a,/J)(*) = J ]im / [Si^\x) TUx, -) + TO^\x) SUx, -) + 

J- «/ — oo ' t ' t 

+TUx, 4) + '^-<^(^' 4) ^n^^(^) ] *(^) (2-7) 

Of course this definition must be understood in the weak (Hilbert) sense. 
Moreover, whenever 5" and T commute (again in the weak (Hilbert) sense), 
the above definition returns the original one. 

3 Multiplying More Distributions 

Up to now we have focused our interest to the multiplication of two distri- 
butions and its possible definitions. This is not enough in many physical 
situations, like, for instance, in the computation of the four-points Green's 
functions in a scalar X(p^ theory, [10]. In this perspective we will now analyze 
possible extensions of the definition (2.4) when more than two distributions 
are considered. In particular we will suggest two different, inequivalent, ap- 
proaches, and we will discuss some examples. Which method has to be chosen 
only depends on which one gives theoretical results in (a better) agreement 
with the experiments (or with the common sense). We will return on this 
point with an example at the end of the next Section. 

In this paper we will consider only commuting distributions. This is an 
useful condition to simplify all formulas. 

The first method we are going to discuss is, in our opinion, the most 
natural one since it does not need any new ingredient for its definition. We 
start with two distributions 5*1 and 5*2. Their multiplication, if it exists, is 
defined by (2.4). Let us now suppose to be interested in defining the product 
of three distributions 5*1, 5*2 and S3 in V'. It is quite natural to consider the 
following quantity 
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liiSi ® S2Yn'''^ S3 + {S, ® Sst'^^ S2 + (^2 ® Sat'^^ Si], (3.1) 



which is certainly well defined for any fixed n, since any term above is the 
product of a C°° function for a distribution. As usual, what may or may not 
exist is the hmit for n ^ 00 of {Si ® ,^2 ® S^)^^''^^ (^f), for any * e V{R). If 
this limit exists we say that the distributions can be multiplied and we put 



It is useful to notice that formula (3.1) would look rather more compli- 
cated without the working hypothesis of the commutativity of the distribu- 
tions. 

Let us now try to define a multiplication between four distribution. In 
this case, of course, we cannot repeat the same steps leading to equation 
(3.2), since the quantity {Si (8) -S'2 <8) 'S'a)^^^) 5*4 would necessarily contain the 
product of two un-regularized distributions. We have to define this multi- 
plication in a different way. This problem can be easily overcome simply by 
coupling the distributions in all the possible ways and then using twice the 
regularization. This implies that the product of four distributions should 
depend on four indices, two a's and two /3's. Explicitly we have: 



whenever this limit exists. To be more explicit, for instance the first term of 
this formula reads 



(^1 ® -S2 ® 53)(.,^) (*) ^ Jim {Si ® ^2 ® ^3)^'^^ 



(3.2) 



{Si ® ^2 ® ^3 ® ^4)(„„„„,,^,) (^) = Jim ^{{Si ® S2t-M(^s, ® S.t^'P^^ + 
+{Si (8) Ss)l^''f''\S2 (8) ^4)!"^'^^^ + (^1 S^)l^'''''\S2 (8) Ss)l^''f^'^ + 
+{ai,Pi)^{a2,P2)}{^), (3.3) 



{Si (8 S2)';^''^'\S3 (8 S^Y^^'^'A (*) 




(3.4) 
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The multiplication of five distributions is now naturally defined in analogy 
with the one in (3.1) and (3.2). We put 

{Si (g) 5'2 O 5'4 (8) S5)^ai,a2,f3uf32)S3 + {Si ^'a O ^'4 (g) S5) (ai,a2,Pi,P2)S2 + 
{S2 ® S's (g) 5*4 ® 5'5)(ai,a2,/3i,/32)'S'l](^), 

whenever the right hand side exists for any '^{x) e T>{R). 

It is clear now how this procedure can be generalized to the multiplication 
of an arbitrary number N of distributions: 

whenever is even we have to proceed like in (3.3), that is we consider all 
the different N/2 pairs of distributions, regularize each pair, and then try to 
remove the regularization. If N is odd, we simply have to multiply one un- 
regularized distribution with the regularization of the even N — 1 remaining 
ones. 

In all the examples discussed in this work we will stick to the situation in 
which all the distributions coincide. In this case all the formulas are strongly 
simplified. Whenever the limits below exist we have: 

{S^S)ia,p) = limSf{x)Sred{x,^) (3.5) 

{S(g)S® S)^a,0) = S Jim (5 ® SYn^f"^ (3.6) 

{S^S^S^ '5)(ai,a„/3„/j.) = „lim (5 St^'f'^^S 5)(r-^^) (3.7) 

{S^S^S^S® S)^a,,a2,PuM = S \im{S S)l^''^'\S S)l^^'^^\ (3.8) 

and so on. The generalization to a bigger number of distributions is straight- 
forward. All the formulas above arc obviously thought in their week forms: 
they must be applied to a generic function ^' e D, like in equation (3.4). 

We now discuss a different proposal which again extends the multiplica- 
tion introduced in (2.4) for two distributions. 
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First of all, let us introduce two complex quantities 01,02, with Oi + 
02 = 2. We modify the original definition (2.4) of the multiplication of two 
distributions by saying that two distributions Si and S2 are A-multipliable if 
there exists a choice of ai and 02, with 01+02 = 2, such that the following 
limit exists: 



4 Loo '^^'^iS^^^ S2,red{x, ^) + 02^1,^6^(2^, ^) Si'^^{x) 



*(x)(ix, (3.9) 



for any * e V{R). 

Of course definition (2.4) turns out to be simply a special case of this 
one when we take oi = 02 = 1. It is interesting to notice that our new 
multiplication depends now not only on a, (3, but also on oi and 02. Of 
course, it may happen that one contribution in (3.9) does not converge for 
n — > 00. In this case, while the multiplication in (2.4) is not defined, the one 
above still exists for a clever choice of Oi and 02. 

The length of the formulas rapidly increases when the number of distri- 
butions to be multiplied grows up. Already for three distributions we need 
to introduce six parameters, an, 012, 013, O21, O22 and 023- whose sum must be 
equal to 6. The A- multiplication of the three distributions is said to exist if 
there exists a choice of the coefficients Ojj's and of the pair {a,P) such that 
the limit below exists: 

{Si 0^2® ^3)f.,;,)(*) = I „lim [oii<)(x)45(x)53,.ed(a;, ^) + 

+oi25i2(x)52,.ed(x, ^)S^^l{x) + aisSiM^, ^)sQ{x)S^^i{x) + 

+0'2lSi^nix)S2,red{x, —)S3^red{x, — ) + a22Si^red{x, —)S'"2)i{x)S2,,red{x, —) + 
TX 71/ Th Th 

1 1 

+a23Sl,red{x, —)S2,red{x, — )4fn(^)] ^i^) dx, (3.10) 

Th Th 

for all e X>. 

In the case of four distributions the number of the coefficients grows up 
to 14, so that it is more and more difficult to correctly keep into account 
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all these contributions. However, the situation drastically simplifies when all 
the distributions coincide. In this case we have symmetry reasons which give 
some extra conditions on the coefficients a. 

For instance, in the case of two equal distributions, from definition (3.9) 
it is evident that we have to take ai — a2 — 1- Hence, this method returns 
the usual result, see equation (3.5). 

Prom (3.10) we deduce that, if Si — S2 — Ss — S, then necessarily 
Oil = ai2 — ai3 =: bi and 021 = 022 = 023 &2, and therefore 61 + 62 = 2. 
Consequently (3.10) becomes now 



{S ® s ® s)f^^p^m 



1 f°° r 1 1 



*(x) dx.(3.11) 



Finally, without going into details, it is possible to prove that for the mul- 
tiplication of four equal distributions we need to introduce three parameters 
ci, C2 and C3, such that 2ci + 3c2 + 2c3 = 7. The multiplication turns out to 
be 

1 1 

{S®S®S® 5)f«,^)(*) = - ]im / [2ci{SlP\x)fSrea{x, -) + 

+?,c,{S^^\x)f{SrM, \)f + 2c,S^i\x){SUx, \)f] ^{x) dx. (3.12) 

The same procedure can be repeated even for a bigger number of distributions 
but we will omit this generalization here since the difficulty grows up very 
fast with the number of distributions. 

Just a comment before ending this Section: in our opinion, this last 
method appears to be less natural than the first one. Nevertheless, we will 
show in Section 4 that it works well in some examples, and its extra degrees 
of freedom may, in turn, be useful in future applications. The main difference 
within the two methods proposed in this Section is that in the first one we 
increase the number of indices a and /3, while in the second one we keep 
this number unchanged but we introduce new extra parameters which were 
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not originally present in the definition we gave in [4]. As we have already 
observed, the preference must be given to that method whose results are 
closer to the experimental data, or the the common wisdom. 



4 Examples: delta functions 

We devote this Section to show how the multiplications defined previously 
work explicitly. In particular, we will show that both methods proposed 
allow to define the product of an arbitrary number of delta function in one 
dimension localized at the same point. The technique we are going to use is 
very much the same as the one used in [4] where two (derivatives of) delta 
functions have been shown to be multipliable. In particular we will need very 
often the well known Lebesgue dominated convergence theorem (LDCT), see 
[11] for example. 

Before starting with the computation of the multiplications we remind 
the readers the expressions of the two regularizations of the delta function, 
[4]. We have: 

5i^\x) = (3 > (4.1) 

and 

Sredix. — ) = a> 0. (4.2) 

Wc begin with considering the first method proposed, eqs. (3.5)-(3.8), 
taking S = 5. We fix first the form of the function generating the delta 
sequence. In this Section we will always assume that ^{x) is the one given 
in (2.6), where F is a given normalization constant (of course wi-dcpcnding) 
and m is an integer which must be taken even so to prevent j\ $(a;) dx to 
be zero. 

The result for [5 ® (^)(Q„fl) is already contained in Section 2, see (2.5). 
Changing a little bit the notation for future convenience, we have, for any 
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* e v{R), 



(5®5)(a,/3)W 




q; > 2/3, 



(4.3) 



where we have defined 




1 ($(t)) 



dt. 



(4.4) 



1 P 



Of course, due to the presence of Ai^2 in (^ ® 5)(q,/3)(^)j we need to take 
m > 2. Otherwise the integral defining Ai_2 would be divergent. 

This result coincides for both the methods proposed: this is obvious since 
the different multiplications introduced in the last Section both generalize 
the multiplication discussed in [4] and refined in Section 2. 

It is very easy to compute the product of three delta functions using our 
recipe; equation (3.6) becomes now 



since $(0) =0 for any m > 0, for any choice of a and f3 in R^. 

Let us now move to the multiplication of four delta functions. The situ- 
ation is no longer so easy. Using (3.7) we have 



Here we are interested to show that there exists a choice of m, ttj and /3j 
for which the limit of the right hand side of this equation exists finite. We 
will show that such a result can be obtained already if we take ai = 02 ='■ 
a and f3i = (32 ='■ P, with some extra conditions on a and (3. We call 
(5 (g) 5 O (5 (g) (5)(„,/3)(^') = {6 6 ® 6 ® 5)(ai,a2A,/32)(*)- Introducing the 
variable t = xn^ in the integral, and using the fact that $(t) has support in 
[— 1, 1], we obtain 



/°° / \ 1 

^ S{x) Si^\x) Sred{x, —)'^{X) dx 

= -$(0) lim n"+'^^(0) = 



(5(8)5«)5(8)5)(,,,„„;3iA)(*) 
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where 

■In U 7r2n2a-5^ (^2 + 1 /^2(a-/3) ) 2 " 

At this point we use the LDCT. In fact, for any a and f3 with 2a > 5(3 we 
find that < 9{t), where g{t) = ^Itp^™"^). Here we have used the 

same notation introduced in [4], and we have called M = sup^g]_ij[ expj^^} 
and L = supjg]_;^ |\E'(t)|. Of course, g(t) is integrable in [—1,1] whenever 
m assumes values bigger or equal to 2. Moreover, the function f^'^\t) 
converges pointwise, whenever 2a > 5/3, to a function /'•"'^•'(t) which is 
equal to zero if 2a > 5(3 and to ^^^^^2^4'^°'* if 2q; = 5(3. In these conditions the 
LDCT can be applied and we get 

(*»*»*«^W,(*)=(„*^^'^^W' f = ^^ (4.5) 

ID, za > op, 

where, of course, m > 2. 

The (®)(a,/3) multiplication of five delta functions is again computed very 
simply. We have 

{S^S^S^S^ <^)(ai,a2,/3i,/32)(*) = 

/oo 11 
6{X) <5f 6red{x, Sredix, —) ^ [x) dx = 0, 

-oo n 

using again the fact that $(0) = whenever m > 0. 

We are now ready to generalize these results: let I be a natural number. 
Therefore, for any * e 'D{R), 

{6^_^^)^a^,..,a„/3,,..A)i'^) = (4.6) 
21+1 

for any choice of and (3i and for $ given by (2.6) with m > 0. On 
the other hand the multiplication of an even number, 21, of delta functions 
may give a non zero (and finite!) result. It depends, in general, on ai and 
(3i with i = 1,2,.., I, see equations (3.5) and (3.7). As we have already 
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discussed for Z = 2, it is actually enough to put ai — a2 — ■■ — ai —'■ a and 
/3i = = .. = Pi =: p. We obtain 



21 

i ^Ai2idm, la = (3l-l)(3 

Obviously, 74;^2i < oo only if m > 2. 

We now move to the second definition of the multiplication we have in- 
troduced in the last Section. We show that also this method gives non trivial 
results. 

We know already that the multiplication of two delta functions is certainly 
well defined, since it coincides with the multiplication obtained following the 
first procedure. In other words, we have 

(S » *)f,«(*) = » = { g'-^''^'"*'' 1 1 1^^^ (4.8) 

and m must be bigger or equal to 2. 

When we consider three delta functions we obtain, from (3.11), 



(5®5®5)f„,^)(*) = 



ijim [6i(5f (a:))25.e.(x,^) + Mi^H^)(5.e<i(:r,^))^ 



*(x) dx.(4.9) 



We will not give here all the details of this computation, which are very 
similar to those discussed above. The steps are, more or less, the same: 
we change the variable in the integrals putting t — xn^, we restrict the 
integration range due to the compact support of ^{t), and then we use the 
LDCT which can be applied under certain conditions on a, (5 and m. For 
example, the first contribution in (4.9) converges to a finite quantity whenever 
a > 3/3 and for m > 1. On the contrary, the second contribution is surely 
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convergent for a > 2/3 and for m > 4. Collecting these results we obtain 
that, for all m > 4, then 



which is, in general, different from the analogous result, {S<S>S<S>S)(^a,i3){'^) — 0, 
obtained using the first method. 

To multiply four delta functions we refer to equation (3.12). For any 
^ e T>{R) we have 



Now we need to estimate, using the usual techniques, three different contri- 
butions: the first is convergent whenever a > 4/3 and for any natural m. The 
second one converges whenever 2a > 5/3 and m > 2. The last term, finally, 
converges if o; > 2/3 and m > 6. We conclude that, for any m > 6, 



It may be worthwhile to notice that the condition on m does not follow 
from the requirement of ^3,2, to be finite. In fact A3 2 < 00 for any natural 
m. It follows from the analogous requirement for Ai^q, which appears in the 
computation of the last contribution in (4.11), the one proportional to C3. 

Of course an extra degree of freedom is present now: the coefficients bi in 
(4.10) and Ci in (4.12) must satisfy only the very weak constraints: 61 + 62 = 2 
and 2ci + 3c2 + 2c3 = 7. But, since 62, C2 and C3 do not appear at all, any 
choice of 61 and Ci is allowed. 

We now generalize the above results. In general we get 




(4.10) 



1 /""^ 1 

{S^S^S^ <^)fa,/3)(*) = 7 Jim / [2c,{5i^\x))%Ux, -) + 

+3c2{5i^\x)n5U^, 4))' + 2cs5i^\x)(5rea{x, ^))^] *(x) dx. (4.11) 




(4.12) 




(4.13) 
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where d is a positive constant and m > 2{l — 1). Again, this constraint 
on m follows from a term which, under these hypotheses on m, a and 
P does not contribute to the final result, that is, the one proportional to 

Remarks .- (a) It is interesting to observe that, for any odd integer n, the 
(®)(a/3) multiplication of n delta functions may be different from zero while, 
the analogous computation made using (®)(a,/3) returns necessarily zero. 

(b) It is straightforward to generalize all the results obtained in this Sec- 
tion even to the multiplication of the derivatives of the delta function. The 
technique is, more or less, the same. We refer to [4] for the details on the 
regularization procedures of the distributions 5^"^ (x) . 

As in [4] we can apply these results to one dimensional physical models 
which describe media with impurities localized in certain fixed points, or to 
the discussion of the classical limit of a certain quantum mechanical situa- 
tion. Let us consider, for instance, a three-particles system described by a 
factorazible wave function 

^'{xi,X2, X3, t) = ^{{Xi, t)mX2, t)mX3, t) 

where 

We know that Pe(a;i, 2:2, Xa) = \^l{xi,Q)\'^\^2{x2:0)\'^\^l{xs, 0)\'^dxi dx2 dx^ 
is the probability of finding at t — particle i between Xi and Xi + dxi, 
i = 1, 2, 3, [12]. In the limit e ^ we get |$-(a;, 0)p 6{x) (for instance in 
V), so that Pe{xi, X2, X3) — > S{xi) S{x2) S{x3) dxi dx2 dx3. Because of this we 
say that e — > corresponds to the classical hmit of the system: in fact each 
particle is sharply centered in a single point. 

We may look, therefore, for the probabihty of finding the three particles 
at the same point x, in this classical limit. Of course simple physical consid- 
erations require this probability to be zero. Therefore, since this probability 
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should be proportional to {5{x))^, we conclude that the natural regulariza- 
tion is the one in (4.6) with I — 1 and for any choice of a and (3, or the one 
in (4.10) with a > 3f3 and m > 4. 

5 Another example: Klein-Gordon model in 
1 + 1 dimensions 

The model of free bosons which we are going to discuss in this Section is 
defined by the following second order differential equation 

and by the equal time canonical commutation relations 

[ip{x,t),^{x',t)] = 

[<f{x,t),^{x',t)] = 0, (5.2) 
[(p{x, t), (p{x', t)] — iS{x — x'). 

Following the notation and the main steps of [13], we expand the solution of 
the Klein-Gordon equation in plane waves, 

\a{k)e'''^-''"^* + at (A;) 6-^'=^+^'"'=* , (5.3) 

where uj^ — ^/k^ + w? and the operators a{k) and its hermitean conjugate 
a){k) are the coefficients of the expansion. They satisfy these canonical 
commutation relations: 

[a{k), a{k')] = [a\k), a\k')] = 0, [a{k), a\k')] = 5{k - k'). (5.4) 

Let us call the ground state of the theory, [13]. This is defined by 
requiring that a{k)^Q — ^k e R. 
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Interesting quantities to compute are the expectation values in of the 
field (p{x,t) and of the product of the field, (p{x,t)(p{x' ,t'). As in the four- 
dimensional situation, even in this simpler model problems arise when we try 
to compute the matrix element of the product (f{x,t)(fi{x,t). In particular 
we observe that 



^-^k\r.-r,)^ (5.5) 

-00 '±T\ijJh 



evident that, in the limit Ty, A^{rx—ry) diverges logaritmically. (Recall 
that in four dimensions the analogous divergence is quadratic.) 

Now we arc ready to discuss the application of the regularizations pro- 
posed to the Klein-Gordon field. In particular we will discuss first the reg- 
ularization of Lp{x,t) when t is considered an extra parameter. This choice 
is necessary, at this stage of knowledge, since the analytical regularization 
has been introduced only in R, while the sequential completion method is 
formulated in i?". The generalization of the analytic regularization to n > 1 
is discussed in [7]. Even if it is easily seen that both the regularization pro- 
cedures work well as far as the smearing of the field is concerned, we will 
also conclude that the multiplication discussed in Section 2 does not allow to 
control the divergence of A_|_(0), even if the time is considered properly and 
not as a parameter. We hope to be able to reconsider positively this problem 
in a future paper. 

We start considering the analytic regularization of the field cp. Using 
definition (2.2) and considering i as a parameter we get for any e > 0, 



while, a straightforward calculation shows that 



where r. 



X 



(xo, —x) and k- {tx — Vy) = uJk{xo — i/q) — k{x — y). It is therefore 




■0 



dk 



Jo 



°° dk 



a{k)e' 



,ikx—iwkt—ke 



a 



y/^nUJk 



— oo 
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where some easy applications of the integration in the complex domain has 
been used. Analogously we get 



^ ' J-oo \/AnUJk Jo \/4:TTUJk 

Therefore the regularized function, cpj.gg{^x, 
can be written as 



\a\k) 6-^*^^+^'^'=* + a{k) 6^'=^-^'^'=* P,{k), (5.6) 

-oo \/A'KUJk '- 



where we have introduced the (even) function Pe{k) = e~'^^9{k) + e'^^9{—k). 
Prom this equation and from (5.3) it is easy to understand heuristically why 
(freg{x,e;t) is called a 'regularization' of </?: it appears evident, in fact, that 
when e — > then (preg converges in some sense to (p. This follows from the 
fact that, when e — > 0, hence Pe{k) — > 9{k) + 9{—k). Therefore, in this hmit, 
this function behaves like the unit function whenever considered 'inside an 
integral'. More precisely, if f{k) is an integrable function, then we have 



oo 



/oo 
dkf(k). 
-oo 



Let us now make this heuristical argument rigorous, showing that (fregi^, ^] t) 
converges to ip{x, t) in the topology of S'{R) whenever e is sent to 0. In par- 
ticular, we are going to show that, in the limit e ^ 0, the following quantity 

/oo 
[iPredix, e; t) - ip{x, t)] Cix) dx ^2), (5.7) 
-oo 

goes to zero. Here ^^i, ^2 are vectors of the Hilbert space, and C,{x) is a func- 
tion in S{R). Using equations (5.6) and (5.3), exphciting the form of Pe{k) 
and introducing the functions ai2{k) = (^1, a{k)'i>2), ciui^) = (^i; ct^(^)*2), 
and the Fourier transform of C(^)) 

1 

ZTT J-oo 
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we can write S^{(p) as the sum of four contributions all, more or less, of the 
same kind. The first contribution, for instance, is proportional to 

Since C,{k) is a function of <S and 012 (fc) is surely well behaved, we can use 
LDCT to conclude that the above integral converges to zero when e goes to 
0. We arrive to similar conclusions also for the other three contributions in 
5e{(p). This implies that ipred{xie]t) converges to (p{x,t) in S' . 

We discuss now the way in which a delta family can be used in the 
regularization of the scalar field. As for the analytic method we consider 
the time as a parameter. Therefore we have 

/oo poo n 

ji^\yMx-y,t)dy^ j_^^q)^{x - ^,t) dq, (5.8) 

where, as usual, we indicate with the function generating the 5-sequence. 
We now prove explicitly that if $ satisfies the following three conditions, then 
ipi^\x,t) ^ ^{x,t) in S': 

ii) ^{x) = $(-a;); 

iii) $ G S{R). 

Incidentally, we observe that such a function generates a delta family of 
Dirichelet type by means of the procedure discussed in Section 2. Condition 
ii), which is not required in the original definition of the functions of this 
family, is only an useful technical requirement. 

Since t is considered as an extra parameter, we need to prove explicitly 
the convergence of ip^J^\x, t) to Lp{x, t). For this reason, similarly to what we 
have done in (5.7), we compute the following limit 

POD 

lim = lim(^i, / [iplP\x,t)-ip{x,t)]C{x)dx^2), 

where ^1,^2 and C,{x) are the same as in 5e (</?). Using the parity of the 
function $ and introducing again the functions ai2{k) and 012 (^), we can 
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write Sn{(p) as the sum of two contributions with the same structure. In 
particular the first term of 5n{(p), Sn,i{'fi), is 

Again, we make use of the LDCT. The procedure is now a bit tricky. First of 
all, since ( belongs to S, as well as it is clear that, Vn e N, the function 

belongs to C^{R). In fact we can write \fn{k)\ < g{k), with 
g{k) = |ai2(A;)C(-^)|(27rM+ 1). 

Here M is the supremum of the function Obviously, since g{k) e C^{R), 
then also /„(A;) e C^{R). 

This implies that, for any e > 0, it is possible to choose a positive quantity 
independent on n, such that i\k\>R^ \fn{k) \ dk < e. Due to the hypothesis 
i) of normalization of the function $, which can also be written in terms of 
its Fourier transform as 27r$(0) = 1, we deduce that, as far as \k\ < R^, fnik) 
surely converges almost everywhere to the function zero. This means that, 
using LDCT 

lim / \f^{k)\dk= f lim \f^{k)\dk^O, 

which also implies that, given e, it exists Ue & N such that, for all n > rie, 
I\k\<Re dk < €. We can conclude that for all e > 0, it exists a natural 

Tie such that, for all n bigger than Ue, 



fn{k) dk 



< 2e. 



An analogous estimate can be performed also for the second contribution of 
5„ ((/?), Sn,2{v)- We conclude that (/^^^^ converges to (p in S'. 
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Defining the following set of functions 

Z = e : y^°° 1, = (5.9) 

we can summarize the above results in the: 

Proposition 2. For all $ e Z the function (/?(f) = $(g) ip{x-q/n^, t) dq 
converges to (p{x, t) in the topology of S'. 

More results on this convergence will be discussed in the Appendix. 

Once we have shown how the regularizations work for the quantum free 
field, we may think to use ip!^\x, t) and cpredi^, t) to eliminate (some) diver- 
gences appearing in the quantum model. For instance we may think that the 
regularization of A+^r^ — Vy) can be made finite for — Ty. Unfortunately 
this is not so. In fact, let us define, as it is natural, 

[A+(r, - r,)](«,^) = i^o, i^ir,) ® 9'(r,))(a,/3)^o), (5.10) 

and let us focuse our attention in particular to [A_|_(0)](a^^). 
We start computing, see (3.5), 

Inif) = (*0, t) ifregix, t)^o), 

and then we discuss the limit of In{f) for n diverging. Considering only 
the non vanishing contributions (a(/c)^'o = '^oa'^ik) = 0), and using the 
commutation relations of the bosonic operators a{k) and a^(/c), we get 

7^((^) = -1 /~ dw^w) r ^e-^'^-Z^'e-^W*^/"", (5.11) 

4:71 J-OO J-OO UJk 

where 6{k) is a function which is equal to 1 for k > and to —1 otherwise. 

Using the fact that, since $(a;) is taken in Z then ^(x) is an even function, 
as well as its Fourier transform, we have 
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Let us introduce now this new set of functions: 



Zo = {^{k) e V{R) : 27r^(0) = 1, ^(A;) = ^(-A;)} . (5.13) 

It is obvious that the Fourier anti-transform, FT"^, of any function in Zq 
belongs to Z, since T> C S. It may be useful to take $ such that l» e Zq, 
since in this way Inif) can be computed easily using numerical techniques. 
With the change of variable q — k/n^, calling again M the supremum of the 
function ^{k), and assuming that the support of is the interval [—1, 1], 
we deduce that 

IM = [' -jJL= e--^-" < M -jJL= e--^-". (5.14) 



In order to get analytic informations on the asymptotic behavior of Ini'p) 
we begin with an easy estimate which shows that the above integral cannot 
be convergent for n ^ oo whenever a > f3. This follows from the following 
analytic estimate: since for q e [0, 1] e~'^"''^ ° > e~^'^ °, it follows that 

dq ^-qnl^-" ^ ^-nP-<^ ^q 



e-«" > e 



, /„2 _|_ rr^ Jo . /„2 i rr^ 



e-"^ log 



m 



Of course, whenever a > j3 the right hand side diverges. This does not 
really imply that also /„(¥?) diverges, as it is clear. Nevertheless it is a very 
strong indication which, moreover, it is also supplemented by the following 
remark: when n — > oo, the first integral in (5.14) behaves, when a > P, 
like Jq ^ ^{q): which can be finite only if l>(g) goes to zero when q ^ 0. 
This is not what we have since $ belongs to Z, so that its value in A; = is 
^. These result suggests that for In{<p) to be converging, P must be chosen 
bigger than a. But also in this case it is not easy to find an analytic estimate 
for the integral in (5.14) proving that lim„_»oo -^n(</7) < oo. For this reason 
we have used numerical procedures to compute this integral, for different 
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choices of a and (3. Unluckily these numerical results seem to show again 
that lim„^oo Iniv) — oo, even if the divergence is very slow. 

Before ending this Section, we briefly comment on the complete rcgular- 
ization of the field, that is the one in which we consider properly t as the time 
coordinate of the field. First of all we notice that, in a certain sense, only 
in this case we are allowed to speak of a canonical regularitazion of the field 
since the general theory says that the two-dimensional convolution ip * 5'£^ is 
a C°° function and that the two-dimensional analytic regularization is an an- 
alj^ic function. The computation of In{^) does not present many differences 
with respect to the situation discussed above and, by the way, the conclusion 
is still the same: we get, for any choice of a and /3, lim^^oo /„(y9) = oo. For 
this reason we believe it is not worthwhile to give here the details of this 
procedure, which are much heavier than those discussed above and, again, 
do not lead to a positive conclusion. 

6 Conclusions 

In this paper we have discussed different generalizations of the multiplica- 
tion of distributions first introduced in [4]. In particular, we have proposed 
possible modifications of the sequential completion method which may be of 
some utility depending on the distribution to be regularized. 

Furthermore, we have introduced two different definitions of multiplica- 
tions oi N > 2 distributions, both of which generalize the definition given in 
Section 2 lor N — 2. Of course, many other generahzations are also possible. 
We have shown how both these definitions can be used to define the multi- 
plications of an arbitrary number of delta functions localized all in the same 
point. A quantum mechanical physical example has been also sketched. 

Finally we have discussed a naive possibility of using our strategy in QFT. 
We have shown that it is possible to regularize the quantum field in many 
ways, but unfortunately we have also shown that the definition proposed in 
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Section 2 does not allow to cancel out the divergence appearing already for a 
free theory. Our future project are therefore to look for some refinement of the 
procedure which allows to overcome this last problem. If this new technique 
can be found, we can also try to extend the theory to four-dimensional models 
and to discuss the divergences coming from the Feynman graphs. The final 
aim is to consider a non abelian gauge theory like QCD, [14]. 

In this analysis we expect that a crucial role will be played by the function 
$ and by the parameters a and (5 which fix the multiplication. They should 
have the same role as the free parameters in renormalization theory, whose 
values are fixed by the experimental data. 



Appendix : A Convergence Remark 

In this Appendix we prove in a different (and easier) way that, whenever 
e Z, then <p^^\x,t) (p{x,t) in S'. 
Let ({x) e S{R). After some easy computation we deduce that 

<fip\x,t)ax)dx - / ^{x,t)ax)dx^ 

-oo J —oo 

_^ {y) - 6{y)) ip{x - t) ({x) dx dy. 

Of course the integral rjly) = /f^ (p{x — y, t) (^{x) dx is continuous in y. 
Using the results in [8], we can conclude that — > for n — > oo. In fact, 
in particular, if is taken positive, then it generates a delta family of 
positive type, so that for any function f{x) continuous in the origin we have 

/oo 
Si^\x)f{x)dx^f{0). 
-oo 

If $ is not positive the same conclusion still holds since ri{y) is also differen- 
tiable in y = 0, see [8]. 
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